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Abstract. Theory of Newtonian dynamical systems admitting normal shift of hy- 
persurfaces was first developed for the case of Riemannian manifolds. Recently it 
was generalized for manifolds geometric equipment of which is given by some regu- 
lar Lagrangian or, equivalently, by some regular Hamiltonian dynamical system. In 
present paper we consider further generalization of this theory for the case, when 
geometry of manifold is given by generalized Legendre transformation. 



1. What is normal shift? Brief historical overview. 

Phenomenon of normal shift is very simple by its nature. Let's consider it 
in three-dimensional Euclidean space M^. Suppose that a is some smooth ori- 
entable surface in M^. At each point p E a one can draw unit normal vector n 

such that n = n(p) would be a smooth 
vector- valued function on a. Let's move 
each point p of cr in the direction of vec- 
tor n(p) to the distance t which is the 
same for all points p G a. Then moved 
points pt would form another surface at as 
shown on Fig. 1.1. Changing parameter t 
we would obtain one-parametric family of 
surfaces. This construction is known as 
Bonnet transformation. 

In Bonnet construction initial surface 
a is transformed by moving each point of 
a. Trajectories of motion in this case are 
straight lines directed along normal vec- 
tors and points of a move along them with a constant speed |v| = 1. Therefore 
parameter t, which is the distance of displacement, can also be interpreted as time 
variable. Bonnet noted that all surfaces at in his construction are perpendicular to 
the trajectories of moving points. For this reason his construction is also known as 
normal displacement or normal shift. 

Basic observation by Bonnet, i. e. orthogonality of surfaces at and shift trajecto- 
ries, gave an impetus for generalization of his construction. This was done by me 
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and my student A. Yu. Boldin in preprint [1] (see also [2], [3]). We replaced straight 
lines in Bonnet construction by trajectories of Newtonian dynamical system: 

r = v, v = F(r,v). (1.1) 

This dynamical system describes the motion of a particle of unit mass according to 
Newton's second law. Vector F in right hand side of (1.1) is a vector of force acting 
on this particle. In order to initialize a shift of surface a we applied initial data 

r| =r(p), v| =v{p)-xi{p) (1.2) 

lt=o lt=o 

to differential equations (1.1). For classical Bonnet transformation |v[ = 1. In our 
construction modulus of velocity vector is not constant. It's initial value on a is 
determined by some scalar function v = y{p). 

Initial data (1.2) with parameter p £ a determine a family of trajectories of 
Newtonian dynamical system (1.1). Points of <j moving along these trajectories 
form one-parametric family of surfaces at ■ Thus we have generalization of Bonnet 
construction. This is the shift of a along trajectories of Newtonian dynamical 
system (1.1). It's clear that this is normal shift for initial instant of time t = Q. 
However, in general, orthogonality of at and shift trajectories gets broken at any 
other instant of time t ^ Q. Only for special Newtonian dynamical systems, i. e. for 
special force fields F = F(r, v), one can keep orthogonality for f 7^ at the expense 
of proper choice of function v = y{p) in (1.2). 

Definition 1.1. Shift of surface cr C along trajectories of dynamical system 
(1.1) determined by initial data (1.2) is called normal shift if all displaced surfaces 

at are perpendicular to shift trajectories. 

Definition 1.2. We say that Newtonian dynamical system (1.1) satisfies normal- 
ity condition if for sufficiently small part of any surface cr C there is a function 
y = u{p) on it determining normal shift of this part of a. 

Words "sufficiently small part" of a in definition 1.2 mean that for any fixed 
point Pq & a there is some sufficiently small open neighborhood of pq in cr, where 
proper function y = v{p) does exist. Let z/q 7^ be some arbitrary nonzero constant. 
Then we can normalize function y{p) at the fixed point pq: 

v{pa) = T^o- (1-3) 

Definition 1.3. Newtonian dynamical system (1.1) satisfies strong normality 
condition if for any surface cr C M^, for any fixed point po G o", and for any 

constant 1^0 7^ there is some open neighborhood of po ^md some smooth function 
y{p) normalized by the condition (1.3) in this neighborhood such that it determines 
normal shift in the sense of definition 1.1. 

Newtonian dynamical systems satisfying strong normality condition form special 
subclass which appears to be interesting object for study. We call them systems 
admitting normal shift of hypersurfaces. In simpler words, these are systems 
capable to implement normal displacement of any hypersurface a in with any 
predefined value y^ of initial velocity. 
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It's obvious that definitions 1.1, 1.2, and 1.3 can be formulated for higher di- 
mensional Euclidean spaces and for Riemannian manifolds as well. One should only 
replace surfaces by hypersurfaces and replace by M" or by arbitrary smooth man- 
ifold M with Riemannian metric g. This was done in papers [4-6]. Instead of (1.1) 
in Riemannian manifolds we write differential equations 

x'=v\ Vtv' = F\x\... ,x'',v\... ,v''), (1.4) 

where i = 1, . . . , n. Here x^, ... , a;" are coordinates of a point p in some local 
chart of Riemannian manifold M, while v^, ... , are components of velocity 
vector V G Tp{M). Vector F G Tp{M) with components F^, ... , F" determines 
force field of Newtonian dynamical system (1.4). 

In papers [5] and [6] we have shown that strong normality condition applied to 
Newtonian dynamical system (1.4) leads to a system of partial differential equations 
for components of force vector F = F(p, v). This system subdivides into two parts: 

- weak normality equations written for n ^ 2; 

- additional normality equations written for n ^ 3. 

Note that n = 2 is lower limit for the dimension of manifold M . Indeed, for n = 1, 
hypersurfaces are points, concept of normal shift in this case has no meaning. Note 
also that additional normality equations are written for n ^ 3. Therefore n = 2 
is exceptional dimension. Theory of dynamical systems admitting normal shift in 
two-dimensional case n = 2 is rather different from that of multidimensional case 
n ^ 3. This fact is reflected in theses [7] and [8]. 

In two-dimensional case strong normality condition is equivalent to weak nor- 
mality equations for force field F. Weak normality equations in this case can be 
reduced to one nonlinear partial differential equation for one scalar function of four 
variables. This equation cannot be solved explicitly in general. However, one can 
construct some special explicit solutions of it (see paper [9] and thesis [8], where 
numerous examples are given). 

In multidimensional case n ^ 3 strong normality condition is equivalent to com- 
plete system of weak normality equations and additional normality equations for 
components force field F. As appeared, in this case one can find an explicit formula 
for general solution of this rather huge system of PDE's (see papers [10], [11], and 
Chapter VII of thesis [7]). Here is this formula: 

_ h{W) v^_^ Vfciy 2^fe^, -IvpJfe 
'~ ' Iv| U ' 1^1 ■ 

Formula (1.5) contain two arbitrary functions W and h, where h = /i(w) is an arbi- 
trary smooth function of one variable, while W = W{x^, . . . , a;", u) is an arbitrary 
smooth function of n -|- 1 variables with nonzero derivative 

»'. = f^^o. 

OV 



By VkW in formula (1.5) we denote partial derivatives 
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while Vi and v'' are covariant and contravariant components of velocity vector v: 

n 

fe=i 

When substituted into (1.5), last (n + l)-tli argument v of Wy and VkW is replaced 
by modulus of velocity vector: v = |v|. 

Division of complete system of normality equations into two parts is not artificial 
by its nature. As appeared, weak normality equations have their own geometrical 
interpretation. In paper [12] normal blow-up of points in Riemannian manifolds was 
considered. There it was found that force field of Newtonian dynamical systems 
admitting normal blow-up of points should satisfy weak normality equations only. 
This means that in multidimensional case n > 3 they could form larger class of 
dynamical systems than those admitting normal shift of hypersurfaces. Examples 
given in paper [13] show that they actually do form larger class. 

Papers [14] and [15] are devoted to the study of global geometric structures asso- 
ciated with Newtonian dynamical systems admitting normal shift of hypersurfaces 
in multidimensional case n ^ 3. This study is based mainly on explicit formula 
(1.5). In addition, one should note that theory of dynamical systems admitting 
normal shift of hypersurfaces was generalized to the case of Finslerian manifolds 
(see Chapter VIII of thesis [7]). Weak and additional normality equations were 
derived. However, explicit formula like (1.5) for this case is not yet obtained. 

During one year after the conference dedicated to Centenary Anniversary of 
I. G. Petrovsky, May 2001, 1 was looking for applications of the theory constructed. 
Results arc represented by papers [16] and [17]. The idea is very simple. It is 
known that in the limit of short waves all wave propagation phenomena can be 
described in terms of rays and beams (geometrical optics is an example). In this 
limit amplitude of scalar wave is described by asymptotic expansion 

oo 

u=y^.e^^^, A ^00, (1.6) 

known as Debye's ansatz (sec [18]). Here in (1.6) function S — S{x^, . . . ,x") is a 
phase of propagating wave. This function satisfies Hamilton-Jacobi equation 

H{x\... ,a;",ViS,... ,V„5) = 0. (1.7) 

Here ViS, . . . , V„S' are components of momentum covector p = VS: 

^'. = V.5=^. (1.8) 

Hamilton fmiction H = H(x^, . . . . . . ,p„) in (1.7) is determined by wave 

operator describing physical properties of medium, where wave propagation process 
occur. Usually it is polynomial with respect to components of momentum covector 
(1.8). But in more complicated cases it may be non-polynomial as well. 
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Note that Hamilton- Jacobi equation (1.7) is a first order PDE with respect to 
phase function S = S{x^, ... , .r"). Its solution is written in terms of characteristic 
lines. In present case they are given by Hamilton equations 

Physically, characteristic lines or trajectories of Hamiltonian dynamical system (1.9) 
are interpreted as rays or beams in wave propagation process. Wave fronts in this 

process are hypcrsurfaccs, where phase is constant. In other words, they arc level 
hypersurfaces of phase function S = S[x^ , . . . ,a;"). These hypersurfaces represent 
the position of wave at various time instants, so we can say that wave propagates 
moving along trajectories of Hamiltonian dynamical system (1.9). However, time 
variable <, with respect to which ordinary differential equations (1.9) are written, 
is not an actual time in wave propagation process. Actual time is proportional to 
the value of phase itself, therefore we can take t = S. As shown in [17], passing to 
this new time variable, we obtain the following modified Hamilton equations: 

Denominator Q, which is often interpreted as kinetic energy, is given by formula 



n O TT 

^=y,j^Pk. (1.11) 



Now suppose that we have Riemannian manifold M. Then we can define modulus 
of momentum covector p = |p| and can consider Hamilton function of special form 
H = W{x^ , . . . , x", |p|). Modified Hamiltonian dynamical system (1.10) with such 
special Hamilton function appears to be equivalent to Newtonian dynamical system 
(1.4) with force field given by the following formula: 

This is the main result of paper [17]. Comparing (1.5) and (1.12) we conclude that 
part of Newtonian dynamical systems admitting normal shift of hypersurfaces in 
Riemannian geometry can be interpreted as the equations of wave front dynamics 
in physics. And conversely, for some wave propagation phenomena we observe 
something like conservation law: wave front dynamics preserves orthogonality 
of wave fronts and rays. However, unlike mathematical theorems, true laws of 
nature cannot be sharply specific, i. e. applicable to some objects and not applicable 
to others. One should cxpcx;t that the above orthogonality law can be generalized 
for all wave propagation phenomena. 

Now suppose that modified Hamilton equations (1.10) are written for arbitrary 
smooth manifold M . In the absence of metric one should revise the concept of 
orthogonality itself. This was done in paper [19]. Note that modified Hamilton 
equations (1.10) define dynamics in cotangent bundle T*M. Therefore at each 
point of trajectory p = p{t) in M we have momentum covector p = p{t). Therefore 
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if trajectory p = p{t) crosses some hypersurface a and if r is a vector tangent to a, 
then we can consider scalar product of r and momentum covector p: 

n 

^=(Pk)=^r'^Pfe. (1.13) 
fe=i 

We see that scalar product (1.13) does not require any metric. If = for all 
T £ Tp{a), then wc say that trajectory p = p{f.) is perpendicular to hypersurface a. 

Let a be some arbitrary orientable hypersurface in M. In the absence of metric 
we cannot define normal vector for a. However, at each point p of a there is normal 
covector n = n(p) perpendicular to a in the sense of scalar product (1.13). It is 
defined uniquely up to a scalar factor: n ^ a . n. We can choose n = n{p) to 
be smooth covector-valued function on a. Then we can define the following initial 
data for ordinary differential equations (1.10): 

x'\ =x\p), pA =v(j>)-ni(p) (1.14) 

lt=o lt=o 

This defines motion of the points of a along trajectories of dynamical system (1.10). 

Definition 1.4. Shift of hypersurface a d M along trajectories of a dynamical 
system in cotangent bundle T*M determined by initial data (1.14) is called normal 
shift if all displaced hypersurfaces at are perpendicular to shift trajectories in the 
sense of scalar product determined by formula (1.13). 

Definition 1.5. Dynamical system in cotangent bundle T*M satisfies strong nor- 
mality condition if for any hypersurface a C M, for any fixed point po G f, and 
for any constant z/q / there is some open neighborhood of po and some smooth 
function v{p) normalized by the condition (1.3) in this neighborhood such that it 
determines normal shift of a in the sense of definition 1.4. 

In paper [19] it was shown that any modified Hamiltonian dynamical system 
(1.10) with arbitrary Hamilton function H (provided ^ 0) satisfies strong nor- 
mality condition. This means that dynamical system (1.10) form background for 
studying strong normality condition. In Riemannian geometry such background is 
given by geodesic flows (these are Newtonian dynamical systems (1.4) with iden- 
tically zero force field F = 0). In paper [19] I considered the following class of 
dynamical systems in cotangent bundle T*M: 

Under sufficiently non-restrictive conditions for Hamilton function H (see [19]) 
differential equations (1.15) can be transformed to the form similar to (1.4): 

x'=v\ i)' = $'(x\... ,«"). (1.16) 

Therefore we say that (1.15) is relative form of Newtonian dynamical system 
(1.16) as related to Hamiltonian system (1.10). Covector Q with components 
Q\, . . . , Qn in (1.15) plays the same role as force vector F in (1.4). 
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Note that Newtonian dynamical system in relative form (1.15) is a dynamical 
system in cotangent bundle T*M. Therefore definitions 1.4 and 1.5 can be applied 
to it. This was actually done in paper [20]. In that paper it was shown that 
strong normality condition for Newtonian dynamical system (1.15) is equivalent to 
a system of partial differential equations for components of covector Q = Q(f^, p). 
This system of partial differential equations subdivides into two parts: 

- weak normality equations written for n > 2; 

- additional normality equations written for n ^ 3. 

Studying these newly derived normality equations is rather interesting problem. 
However, below we go further and we consider more general situation, when Hamil- 
tonian and/or Lagrangian dynamical system in M is not given. 

2. The idea of further generalization. 

In the absence of special geometric structures (like Ricmannian metric or Hamil- 
ton function) the only way of defining Newtonian dynamics in M is given by the 
equations (1.16). Can we use scalar product (1-13) for to define normal shift in this 
case. The answer is yes, provided we have some way to determine momentum co- 
vector p ! For instance, momentum covector p can be given explicitly as a function 
of dynamic variables p and v. In local chart this looks like 

' Pi =pi(a;\... ,a;",i;\... ,i;"), 

< (2.1) 

. Pn=Pn{x^,--- ,a;",u\... ,t;"). 

If coordinates a;^, . . . , a;" of the point p € M are fixed, then functions (2.1) express 
n variables pi, ... , p„ through another set of n variables v^, ... , u". Therefore 
one can treat (2.1) as coordinate representation of a map 

X:TM^T*M. (2.2) 

This map is an analog of well known Legendre transformation (see [21]), for which 
functions (2.1) arc determined by Lagrange function L = L{p.v). Below wc shall 
call it generalized Legendre transformation and for the sake of convenience 
we shall assume this map (2.2) to be diffeomorphism. 

Using generalized Legendre transformation, wc can transform Newtonian dy- 
namical system (1.16) to new dynamic variables p and p forming point q = {p,p) 
of cotangent bundle TM. As a result we get the following differential equations: 

(x' = V''{x\... ,x",pi,... ,pn), ^^^^ 
\pi =6i(a;\... ,a;",pi,... ,p„). 

Here V^, . . . , are functions that implement inversion of the map (2.2): 

' = V^{x^,... ,a;",pi,... ,pn), 

< (2.4) 

_ v" = V"{x\... ,x^,pu... ,Pn). 
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They define kinematic part of the equations of Newtonian dynamics (2.3). Their 
role is similar to the role of Hamilton function in (1.15). Functions @i, . . . ,0„ 
define dynamical part of the equations (2.3), they play the same role as functions 
Qi, ... , Qn in (1.15) and functions . . . , $" in (1.16). 

Note that Newtonian dynamical system written as (2.3) is a dynamical system 
in cotangent bundle T*M. Definitions 1.4 and 1.5 are applicable to it. Therefore we 
can start new theory of dynamical systems admitting normal shift of hypersurfaces 
for manifolds equipped only with generalized Legendre map. Our main goal here 
is to derive weak and additional normality equations for this case. 

3. Newtonian dynamics in manifolds 
and associated geometric structures. 

Let's consider Newtonian dynamical system in some smooth manifold M. In 
the absence of metric we cannot use covariant derivative in (1.4). Therefore we 
write it as (1.16). This is the system of first order ordinary differential equations 
given by Newtonian vector field $ in tangent bundle TM: 

^ f) ^ f) 

1=1 i=l 

Let q = {p, v) be a point of tangent bundle TM. Then canonical projection 
TT : TM M maps it to the point p G M. Canonical projection induces lin- 
ear map TT* : Tq{TM) — + Tp{M). Applying this map to vector (3.1) at the point 
q G TM, we get velocity vector v at the point p e M: 

7r.(*)=v. (3.2) 

The equality (3.2) can be taken for the definition of Newtonian vector field (3.1). 
Linear map tt* in (3.2) acts from Tq{TM) to Tp(M). However, one can consider a 

map acting in opposite direction. 

Definition 3.1. Suppose that for each point q of tangent bundle TM we have 
linear map /: Tp{M) Tq{TM) from tangent space Tp{M) at the point p = Tr{q) 
to tangent space Tq{TM) at the point q. This construction is called a lift of vectors 

from M to tangent bundle TM . 

Definition 3.2. Lift of vectors / from M to TM is called smooth lift if it maps 
each smooth vector field in M to a smooth vector field in TM . 

Definition 3.3. Lift of vectors / from M to TM is called vertical lift if compo- 
sition TT* o / is identically zero: tt* ° / = 0. 

Definition 3.4. Lift of vectors / from M to TM is called horizontal lift if 

composition tt* » / is the field of identical operators on M, i. e. tt* ° / = id. 

Each smooth manifold M possesses canonical vertical lift of vectors from M to 
TM. It is defined as follows. Let X be some vector field in M and let Xp be its value 
at the point p e M. Then we can define one-parametric group of diffeomorphisms 
ipt in TM that maps point q = {p, v) to the point ipt{q) = (p, v+f -Xp). Vector field 
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Y in tangent bundle TM associated with one-parametric group of diffeomorphisms 
ipt is taken for the result of lifting vector field X: 

Y = w(X). (3.3) 

Let's write (3.3) in local chart. Suppose that X is given by its coordinates: 

X = ^^-^ + ... + ^"-^. (3.4) 

Then result of applying canonical vertical lift w to X is given by formula 

Y = .(X)=X^.^ + ... + X".A. (3.5) 

Looking at (3.4) and (3.5), we see that for each point q of tangent bundle w is 
injective map that maps tangent space Tp{M) onto the vertical subspace 

(TM) =Ker7r*. (3.6) 

Let h be some horizontal lift of vectors from M to TM . For each point q of tangent 
bundle h is also an injective. Let's denote by Hq(TM) its image 

Hq{TM) = lmh. (3.7) 

Subspaces (3.6) and (3.7) are transversal and complementary to each other: 

Tq(TM) = Hg{TM) ® Vg{TM). (3.8) 

Lemma 3.1. Defining horizontal lift of vectors from M to TM is equivalent to 
fixing horizontal subspace Hq{TM) complementary to vertical subspace Vq{TM) in 
Tq{TM) at each point q of tangent bundle TM. 

Suppose that we have some horizontal lift h. Lot's denote by Ei, ... , E„ base 
of coordinate vector fields for some local chart of the manifold M: 

Ei = ^,-,E„ = ^. (3.9) 

Due to (3.8) and due to the equality tt* <> /i = id, applying hto Ej, we get: 

ME.) = ^-Er.^-^- (3.10) 

Quantities = (x^, . . . , x", u^, . . . , w") represent horizontal lift h in local chart. 
If they are smooth functions of their arguments, then h is smooth lift. Under the 
change of local chart in M these quantities are transformed as follows: 

= E E c + E E ^ (3.11) 

m=l a=l m=l r=l 
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Here S^, T°-, and T'"' are components of direct and inverse transition matrices 
(Jacobi matrices) for the change of local coordinates: 



S] = Ti = 1^. (3.12) 



In arbitrary smooth manifold there is no canonical horizontal lift of vectors. 
However, if Newtonian dynamical system (1.16) is given, then there is horizontal 
lift /i, canonically associated with it (see [22] and references therein). Let $ be 
Newtonian vector field (3.1) and let X be some arbitrary vector field in TM. Then 
we associate with X the following vector field denoted by H(X): 

H(X) = X+[^;c7r.(X),^]-^;c7r.([X, ^]) _ ^^^^^^ 

Here by square brackets we denote commutator of two vector fields. 

Lemma 3.2. For any smooth junction (f and for arbitrary smooth vector field X 
in TM the equality H((y9 • X) — ■ H(X) is fulfilled. 

Proof. By direct calculations we find that [i^ • X, = • [X, — ■ X. Here 
by $(y9 we denote derivative of function along vector 

Composition wott* acts as linear operator at each point q G TM. Therefore 

WoTT^{[ip ■ X, *]) = ip • wo7r*([X, *]) - • M;o7r*(X), (3.14) 
wo7r*((/?-X) = (/?-wo7r*(X), 

[w<'n*{(p ■ X), = • [M;o7r*(X), 4>] - • wo7r,(X). (3.15) 

When subtracting (3.14) from (3.15), terms containing ^ip cancel each other. This 
proves required equality H((^ • X) = Lp ■ H(X). □ 

Formula (3.13) defines an operator H acting on vector field X and yielding 
another vector field H(X). Commutators in (3.13) contain differentiation with 
respect to X. Therefore one might expect H to be differential operator. Lemma 3.2 
means that H is non-differential linear operator acting pointwise at each point q of 
tangent bundle TM. 

Lemma 3.3. Linear operator H: Tq{TM) Tg{TM) defined by formula (3.13) 
satisfies the equalities H » H = H, tt, ° H = tt*, and H » w = 0. 

Proof. In order to prove required equalities we shall use direct calculations in local 
coordinates. Let's take an arbitrary vector field X in TM represented as 

i=l i=l 
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Newtonian vector field * is represented by formula (3.1). Hence for [X, we have 



fe=l fe=l i=l ^ ' 

i=l k=l ^ ' i=l k=l ^ ' 

Applying operator wott* to (3.16) and to the above expression, we obtain 

»,.^.(X)=£x'.g-, (3J7) 

» . ,.(1X, *]) = t ■ - 1 1 («' + 1^ ) • (3.18) 

fe=l i=l fe=l ^ ' 

Now, relying upon formula (3.17), we calculate the commutator [■u;o7r»(X), $]: 

[„...(x),*l = E^'-^ + EE^'Sr-^- 

fe=l k=l i=l 

(3.19) 



Then subtract (3.18) from (3.19). As a result we obtain the following equality: 

n o 

[w;o7r*(X), *]-«;o7r4[X, *])=^x'=.-^ + 

fc=i 

Qyt Qyk Qyk 
k=l i=l k=l 

And ultimately, for the result of applying operator H to vector X given by formula 
(3.16) we obtain the following expression: 

fe=l 8=1 i=l 

Looking at the expression in right hand side of (3.20), we see that the required 
equalities HoH = H, tt* oH = tt*, and How = for operator H appear to be 
obvious. Lemma 3.3 is proved. □ 

First equality H o H = H means that H is an operator of projection. It projects 
Tq{TM) onto some subspace Hq{TM), where Hq{TM) = ImH. Each projection 
operator breaks the space into direct sum of its kernel and its image: 

Tg(TM) = Ker He ImH. (3.21) 
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Second equality tt, o H = tt* yields Ker H C Ker tt* , while third equality H » w = 
means that Imw C KerH. Taking into account (3.6), we obtain that kernel of 
operator H coincides with vertical subspace Vq{TM): 

KerH = Ker TT* = Imw = y,(rM). (3.22) 

Due to (3.22) the equality (3.21) can be rewritten as (3.8). This means that operator 
H defines horizontal subspace Hq{TM) = ImH in TqiTM). Applying lemma 3.1, 

we derive the following theorem. 

Theorem 3.1. Each Newtonian dynamical system (1.16) in smooth manifold M 
generates horizontal lift of vectors from M to TM canonically associated with it. 

This result is not new (see [22] and references therein). We reproduced it here 
for the sake of completeness of our consideration. 

Let h be horizontal lift canonically associated with Newtonian dynamical system 
(1.16). Then H = /iott* and h^Ei) = hoTr^{d/dx') = li{d/dx'). Applying formula 
(3.20) to vector field X = Ej, we now obtain the following equality: 

r) " 1 (9(1>'' 8 
ME.) = H(a/a.0 = ^+E2^-^- (3-23) 

Comparing (3.23) with (3.10), for components of h we derive explicit formula: 

= -- (3.24) 

4. Extended tensor fields and covariant differentiations. 

Concept of extended tensor fields is closely related to Newtonian dynamical 
system. Indeed, if we look at right hand side of (1.4), we see that F = F(p, v) is 
a vector- valued function with values in tangent spaces Tp{M). However, it doesn't 
fit standard definition of vector field in M since it depends not only on a point 
p G M, but also upon velocity vector v at this point. Both p and v form a point 
q = {p,v) of tangent bundle TM. Function F = F(p, v) is an example of extended 
vector field. Vector-function V = V(p, p) in (2.3) with components (2.4) is another 
example. Its arguments form a point q = (p, p) if cotangent bundle T*M. 

In order to formulate definition of extended tensor fields in general let's consider 
the following tensor product of tangent spaces and their dual spaced: 



r;(p, M) = Tp(M) (E)...<E) Tp{M) t;(m) ® . . . r;(M) 



Definition 4.1. Extended tensor field X of type (r, s) in v-rcprcsentation is a 
tensor- valued function with argument q = (p, v) in tangent bundle TM and with 
value X(g') in tensor space TJ(p, M), where p = n{q). 
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Definition 4.2. Extended tensor field X of type (r, s) in p-representation is a 
tensor-valued function with argument q — {p, p) in cotangent bundle T*M and 
with value X(g) in tensor space TJ(p, M), where p = n{q). 

At first let's consider extended tensor fields in v-representation. Denote by 
Tg{M) the set of smooth extended tensor fields of type (r, s) and take the sum 

OO OO 

t(m) = 00t;(m) (4.1) 

In T{M) we have all standard tensorial operations like summation, multiplication 
by scalars, tensor product, and contraction. Direct sum (4.1) possesses structure 
of graded algebra over the ring of smooth scalar functions in TM, we denote this 
ring by S^(TM). Algebra T(M) is called extended algebra of tensor fields. 

Definition 4.3. A map D : T(M) T(M) is called differentiation of extended 

algebra of tensor fields, if the following conditions are fulfilled: 

(1) concordance with grading: D{T^{M)) C TJ(Af); 

(2) M-linearity: D{X + Y) = ^(X) + D{Y) and D{XX) = XD{X.) for A G R; 

(3) commutation with contractions: I?(C(X)) — C(D(X.)); 

(4) Leibniz rule: ^(X ® Y) = D{X.) ® Y + X ® D{Y). 

Let's denote by £)(M) the total set of differentiations of extended algebra of ten- 
sor fields T(M). It is easy to see that it possesses the structure of 5^(rM)-module. 
The set of extended vector fields Tq{M) is equipped with the same structure of 
module over the ring ^{TM). This coincidence motivates the following definition. 

Definition 4.4. Covariant differentiation V in the algebra of extended tensor 
fields T(M) is a homomorphism of 5^(TM)-modules V : T^{M) D{M). Image of 
vector field Y under such homomorphism is called covariant differentiation along 
vector field Y. It is denoted by Vy- 

In Chapter III of thesis [7] it was shown that each covariant differentiation V in 
extended algebra of tensor fields is associated with some lift of vectors from M to 
TM. It is called horizontal covariant differentiation (or vertical covariant 
differentiation) if corresponding lift of vectors is horizontal (or vertical). In each 
smooth manifold there is canonical vertical covariant differentiation associated with 
canonical vertical lift of vectors w. We denote it by V. In local chart this covariant 
differentiation is given by the following formula: 

VmX]'-Y = —iiiiiii. (4.2) 

Due to (4.2) this covariant differentiation is also called velocity gradient. 

To define horizontal covariant differentiation, apart from horizontal lift of vectors, 
one need some extended affine connection T. Its components are given by formula 

n 

VE,E,=^r^^.Efe. (4.3) 
fe=i 
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Here Ei, . . . , E„ arc coordinate vector fields (3.9) in M . Unlike traditional affinc 
connection, components of extended affine connection T^^ depend on double set of 
arguments, i. e. on coordinates of point p G M and on components of vector v: 

T^^=r':j{x\... ,x-,v\... ,vn- (4.4) 

Under the change of local chart quantities (4.4) are transformed as follows: 

^io = E E E ^/ f ™ + ^ 5i ^ (4.5) 

m=l a=l c=l m=l 

(matrices S and T are defined in (3.12)). If components of horizontal lift h in (3.10) 

and components of extended affinc connection F in (4.3) are given, then horizontal 
covariant differentiation V in local coordinates is given by the following formula: 



a=l 6=1 

r n s n 

I \ ^ \ ^ -rik -v-ii-.-ak-.-ir _ \ ^ \ ^ pfcfc Y'l--- 

^ mak 3s mjk : 



(4.6) 



■m.3k ji---bk 
fe=l afe=l fe=16fc=l 



Note that F*^ and Ff^ are two independent sets of parameters in (4.6). The only 
condition is that they should obey transformation rules (3.11) and (4.5) respectively. 
However, one can prove the following two lemmas. 

Lemma 4.1. If horizontal lift of vectors from M to TM is given and ifV^ are its 
components, then there is symmetric extended affine connection with components 

1 9F'= 1 9Ff 

Lemma 4.2. // extended affine connection with components T^j is given, then there 
is horizontal lift of vectors from M to TM with components 

n 
3=1 

These two lemmas are proved by direct calculations on the base of formulas 
(3.11) and (4.5). Applying (4.7) to (3.24), we obtain 

, 1 a^^fe 

2dv'dvi' ^ ' 

As a result we have proved the following theorem. 

Theorem 4.1. Each Newtonian dynamical system (1.16) in smooth manifold M 
generates extended affine connection with components (4.9) canonically associated 
with this dynamical system. 

Usually equalities (4.7) and (4.8) cannot be fulfilled simultaneously. In all previ- 
ous papers we defined F^ through F^^- by means of formula (4.8). This is equivalent 
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to the equality ViV^ = 0. In present paper we keep this tradition. Therefore 
formula (4.6) for horizontal covariant derivative looks like 



a-v'ii...ir n n n riY^'^" 

a=l b=l c=l 



/ -J / -J in ilk ji js / V / ^ ^1'' jk 3 



(4.10) 



= 1 afe=l fe=l 6,^=1 



Theory of differentiation for extended tensor fields in p-representation is a little 
bit different. Here we also can define extended algebra of tensor fields by 
means of direct sum (4.1). It is algebra over the ring of smooth functions in cotan- 
gent bundle ^{T*M). Definition 4.3 remains unchanged. Definition 4.4 is replaced 

by the following two definitions. 

Definition 4.5. Covciriant differentiation V in the algebra of extended vector 
fields T(M) is a homomorphism of 5^(T*M)-modules V: T^{M) D{M). Image 
of vector field Y under such homomorphism is called covariant differentiation along 
vector field Y. It is denoted by Vy- 

Definition 4.6. Contravariant differentiation V in the algebra of extended 
vector fields T(M) is a homomorphism of S^(T*M)-modules V : T^{M) £)(M). 
Image of covector field q under such homomorphism is called contravariant differ- 
entiation along covector field q. It is denoted by Vq. 

Instead of velocity gradient (4.2) here in p-represcntation we have canonical 
vertical contravariant differentiation V. It is given by formula 

^j^ii... i,. 
OPm 

Due to (4.11) contravariant differentiation V is called momentum gradient. In- 
stead of canonical vertical lift of vectors in p-representation for each point q = {p, p) 
of cotangent bundle T*M we have injective map w : T*{p) Tq{TM) associated 
with differentiation V. Suppose that covector q is given by its coordinates: 

(i = qi-dx^ + ... + qn- dx"-. (4.12) 

Applying w to q, we obtain a vector Y = w{q) given by the following expression: 

d d 
Y = «;(q) = 91 • — + ... + ?„ •—. (4.13) 

dpi dpn 

Formulas (4.12) and (4.13) are similar to formulas (3.4) and (3.5). However, there 
is invariant way to determine vector Y. It is generated by one-parametric group of 
diffeomorphisms (fit in T*M which is defined as follows. If g = {p, p) is a point of 
cotangent bundle T*M, then (pt{q) — (p, P + * ■ q)- 

Injective map w : T*{p) Tq{TM) defined just above satisfy the equality 
-K^ow = Q. Its image coincides with vertical subspace in Tq(T*M): 



lmw = Vg{T*M) = Key:n^. 
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Each horizontal covariant derivative in p-representation is given by formula 

- + 1. 1. + 

r n s n 
I \ ^ \ ^ pjfc -v^ii... Ofc... i,. _ \ ^ \ ^ pb)! v^il ir- 

fc=l afc = l fc=16(. = l 

Here Fj^ and F^^- are components of some horizontal lift of vectors from M to T*M 
and some extended affine connection respectively. Quantities Ff^ are defined by the 
equality (4.3). They satisfy the same transformation rule (4.5) as quantities (4.4). 
However, they differ from (4.4) since they have another set of arguments: 

T%=V%{x\... ,x^,p^,... ,pn). (4.14) 

Horizontal lift of vectors h from M to T*M satisfies the equality w^oh = id. It's 
components are determined by the following formula: 

K=l 

Under change of local chart quantities Tij are transformed as follows: 

= E E + 5: 5^ (4.15) 

a=l c=l k=l m=l 

Comparing formula (4.15) with (4.5), one can prove two lemmas which are similar 

to lemma 4.1 and lemma 4.2. 

Lemma 4.3. If horizontal lift of vectors from M to T*M is given and ifTij are its 
components, then there is symmetric extended affine connection with components 

pfc = (4 16) 

dpk ^ ^ 

Lemma 4.4. // extended affine connection with components Tfj is given, then there 
is horizontal lift of vectors from M to T*M with components 

n 

fe=i 

Usually equalities (4.16) and (4.17) cannot be fulfilled simultaneously. In previ- 
ous papers we defined F^ by means of formula (4.17). This is equivalent to 

ViPj = 0. 

In present paper we keep this tradition. Connection components F^^- will be im- 
ported to p-representation by means of generalized Legendre transformation (2.2) 
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(see below). Then formula for horizontal covariant differentiation looks like 

a^l b^l c^l 



s n 



(4.18) 



I \ ^ \ ^ pifc y-ii... ak...ir _ \ ^ \ ^ p&fc vii ir 

'^2-^ mar, ^ji js Z-^ 'mjk ji---bk---js' 

k=l ak = l k=lbk—l 



Formulas (4.10) and (4.18) mean that we do not use horizontal lift with components 
(3.24) directly, but only for deriving extended connection with components (4.9). 

5. Generalized Legendre transformation. 

In section 2 we learned that some extended affine connection and some horizontal 
lift of vectors from M to TM are canonically associated with Newtonian dynam- 
ical system (1.16). As for dynamical system (2.3) in T*M, I don't know similar 
direct constructions for T^j and Tij in this case. Therefore now we shall import 
these quantities from v-representation to p-representation by means of generalized 
Legendre transformation A. For in (4.14) we write: 

r^,. = r^,oA-i. (5.1) 

This equality (5.1) means that we simply change arguments in (4.4) by substituting 
functions (2.4) foY v^, ... , v". 

Quantities T^j in v-representation are given by formula (4.9). Using formula 
(5.1), one can find their counterparts F^^- in p-representation. However, it would be 
natural to express them through functions V^, ... , V" and Oi, ... , 0„ that deter- 
mine dynamical system (2.3). Differential equations (2.3) describe the same dynam- 
ics as differential equations (1.16), but in other variables x^, . . . , x", pi, . . . , Pn- 
Therefore we can calculate time derivative of velocity vector v as follows: 

j=l 1=1 

Comparing (5.2) with second part of equations (1.16), we obtain the equality 



*'-^-=f:f?^-+i:^e,. (5.3) 

i=l i=l 



Note that partial derivatives dV'^/dpr in (5.3) form Jacobi matrix for transforma- 
tion given by functions (2.4), when a;^, . . . , a;" are fixed. Let's denote 

ff- = — = V'-V\ (5.4) 

dpr 



Formula (5.4) determines components of an extended tensor field g. Tensor g plays 
the role similar to that of dual metric tensor in Riemannian geometry. It is non- 
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degenerate: detg 7^ 0, but, in general, it is not symmetric: g^^ 7^ g^'^. By gij we 
denote components of inverse matrix for (5.4), i.e. 



T.9''9j'' = ^l T.9iJ9"' = St (5.5) 

They form another extended tensor field, which traditionally is denoted by the same 
symbol g. Though being non-symmetric extended tensor field, it is direct analog 
of metric tensor in Riemannian geometry. 

Let / = f{x^, . . . , x", v^,. . . , t;") be some function in v-representation for some 
local chart of M. This might be component of extended tensor field, component of 
extended connection, or coordinate representation of some other geometric object. 
Then /° is its p-representation. Using exphcit form (2.4) of the map A~^, we 
can derive the following transformation rules for partial derivatives: 

i=l ^ ' 

Taking into account (5.4) and (5.5), we can rewrite (5.6) as 

d! .-1 v„ g(/°A-^) 
fe=i 

Now let's apply (5.7) to the function / = $''(a;^, . . . ,a;",w^, . . . ,u") in (5.3): 

r=l \m=l m=l / 

Applying (5.7) once more and taking into account (4.9) and (5.1), we obtain 

r—1 s—1 \ \m—l m—1 // 

For the sake of convenience this formula should be slightly transformed: 



\m=l m=l 



r=l s=l a=l \m=l rn=l 



(5.^ 



Looking at (5.8), it is obvious that Tfj keep symmetry T^j = Tj^ in p-representation 
as well. It is not surprising due to (5.1), but now it is explicitly evident. 
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Theorem 5.1. Each Newtonian dynamical system (2.3) in cotangent bundle T*M 
generates extended affine connection with components (5.8) canonically associated 
with this dynamical system. 

Theorem 5.1 reformulates previous theorem 4.1 with respect to the same Newto- 
nian dynamics, but transferred from tangent bundle to cotangent bundle by means 
of generalized Legendre transformation (2.2). 

6. Regularity condition. 

Now wc arc ready to study Newtonian dynamics given by the equations (2.3) 
without referring to its v-representation (1.16). If we consider trajectory of dynam- 
ical system as a curve p = p{t) in M, then velocity vector on this trajectory is given 
by the vahic of extended vector field V = V(p, p). Let's define another extended 
vector field W = W(p, p) with components 

n 

W' = ^W'Vpr (6.1) 
r=l 

and consider scalar product of this vector field W and momentum covector p: 

n 

n = {ip\VJ)=Y,PsW\x\... ,x^,v\... ,v^). (6.2) 

s=l 

Extended scalar field f2 = f2(p, p) in (6.2) is somewhat like kinetic energy in me- 
chanics (compare (6.2) with (1.11) and (1.13)). In mechanics both momentum and 
kinetic energy represent the "amount of motion" stored in moving object. This 

motivates the following definition. 

Definition 6.1. Generalized Legendre map (2.2) given by components of extended 
vector field V in (2.4) is called regular if 

(1) it is diffeomorpliic: 

(2) V{p, p) = is equivalent to p = 0; 

(3) f}(p,p)^Oforp^O. 

Third part of regularity condition means that vector W = W(p, p) is transversal 
to null-space of momentum covector p. Therefore one can consider operator-valued 
extended tensor field P with the following components: 

p; = s}-^- (6.3) 

Operator P = P(p, p) with components (6.3) is a projector onto null-space of 
momentum covector p along vector W = W(p, p). It is defined everywhere in 
T*M except for those points q = (p, p), where p = 0. 

7. Force covector. 

Let's consider Newtonian dynamical system (2.3) in smooth manifold M. Sup- 
pose that M is equipped with symmetric extended aflane connection F. This might 
be connection (5.8) canonically associated with dynamical system (2.3) or any other 
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symmetric extended affinc connection which is not related to (2.3) at all. In both 
cases one can introduce extended covector field Q with components 



j=l fe=l 

Then, using Qi, . . . , Qn, one can write differential equations (2.3) as follows: 
(x' = V'{x'^,... ,a:",pi,... ,pn), 

\ ^tPi = Qi{x^,... ,x"-,pi,... ,Pn). 



(7.1) 



(7.2) 



Components of covector Q in (7.2) play the same role as components of vector F 
in (1.4). Therefore covector field Q is called force covector or force field of 
Newtonian dynamical system (7.2). 

8. Weak normality condition and weak normality equations. 

Let's consider one-parametric family of trajectories of Newtonian dynamical sys- 
tem (7.2). This is a family of parametric curves q = q{t, y) in T*M, where t is time 
variable and y is additional parameter. In local chart this family of curves is rep- 
resented by the following set of 2n functions 



' = x^{t,y), 



x^{t,y), 



' v\ =Pi{t,y), 



, Pn =Pn{t,y), 



(8.1) 



Differentiating first part of these functions (8.1) with respect to additional param- 
eter y, we obtain vector-function r = T{t, y) with components 



dy ' 



(8.2) 



Functions pi, ... , p„ in (8.1) arc components of covector-function. Therefore we 
should apply covariant derivative to them: 



^^-v.p^ = f:-EEr^.^^-^- (8-3) 
^ j=i fe=i 

As a result we get covector-function £ = ^{t,y) with components (8.3). Vector r 
is called variation vector or, more exactly, vector of variation of trajectory, 
while covector 4 is called covector of variation of momentum. 

Components of both functions T(t, y) and |(i, y) satisfy a system of ordinary dif- 
ferential equations with respect to time variable. In order to derive these equations 
let's differentiate equations (7.2) with respect to parameter y. This yields 



n n 
fe=l k=l 



(8.4) 
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+ E E E - E E ^tl ps )-r' + 

k=l \j=l 5=1 j=l 5=1 / 

n n n n n \ ' J 

+ E E E ■^k = J2 ^fc^^ ■ + E ^"^^ • 

fe=l j=l s=l fe=l fe=l 

Here in (8.5) we used components of curvature tensors R and D. First is an analog 
of standard curvature tensor of Riemannian geometry, it is given by formula 

m=l m=l 

8.6 



+ E E ■'^'"^ E E -"^"ij 



m=l a=l m=l a=l 

Second is a tensor of dynamic curvature, it is nonzero only for extended connections, 
when V\- do actually depend on components of momentum covector: 

^ = (8-7) 

In the next step we consider scalar product of vector r and momentum covector 
p. This scalar product introduced by formula (1.13) determines so called deviation 
function ip. Like t = T{t,y) and ^ = ${t,y), this is a function of time variable 
t and additional parameter y. In general case deviation function (1.13) satisfies 
linear homogeneous ODE of the order 2n (see theorem 6.1 in [20]). But here we 
consider special case determined by the following definition. 

Definition 8.1. We say that Newtonian dynamical system (2.3) satisfies v^^eak 
normality condition if for each its trajectory q = q(t) there is some second order 
homogeneous linear ordinary diff'erential equation 

^ = Ait) if + Bit) if (8.8) 

such that all deviation functions on the trajectory satisfy this differential equation. 

Saying "all deviation functions", in definition 8.1 we imply that each trajectory 
q = q{t) can be included into one-parametric family of trajectories by various 
possible ways. Each such inclusion defines some variation vector r = r(i) and 
corresponding deviation function (f = (p{t) on that trajectory. Functions A{t) and 
B{t) in (8.8) depend on the trajectory q — q{t), but they do not depend on how 
this trajectory is included into one parametric family of trajectories. 

Let's calculate time derivatives of deviation function (1.13). For first order time 
derivative cp we obtain the following expression: 

n n 

^ = Vt^ = Y, Pi + J2^' VtPi. (8.9) 

i=l i=l 

Then we substitute (8.4) for Vtr' and (7.2) for VtPi into (8.9). This yields more 
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detailed expression for first order time derivative of deviation function: 

n / n \ n / n \ 

fc=i \ i=i / fc=i \ i=i / 

Applying Vt to (8.10), we derive formula for second order time derivative (p: 

n / n \ n / n \ 

A;=l \i=l / k=l \i=l / 

n n / n \ n / n \ 

+ EE E^'-^fe^>' + ^-^'^ y-.r'= + ^ 5^VfeF^Qj •r'= + 

n n / n \ n / n \ 

n n / n \ n n / n \ 

+ EE E^'-^'^'^'O^^'^'^+EE E'^'^^'^^O^^'^'^- 

fe=lr=l\i=l / fe=lr=l\i=l / 

Let's substitute (8.4) for Vtr*^ and (8.5) for Vt^fc into the above equality. It's easy 
to note that resulting expression for will be of the form 

n n 

^ = E"'^fc + E'^fc^'- (8.11) 
fe=i fe=i 

Here a^, ... , a" are components of extended vector field a given by formula 

n n n n n 

a" = Y,^''v'Qi+J2Yl p» + E E V'- v'^y^ ft 0. + 

i=l r=l i=l r=l i=l 

n n n n n n 

+ 1212^''^'Pi^'Qr-J21212Yl ^rj V'^F^ ft ft + (8.12) 
r=l i=l r=l 5=1 i=l j=l 

n n n 

+ E E Vrl^^ft V'^y'^ + J2 V'^'^ Qr. 
r=l i=l r=l 

Similarly, quantities /3i, ... ^ f3n are components of extended covector field (3: 

n n n n 

r— 1 2—1 r— 1 2=1 

71 n n n 

+ E E V" V;tl/* ft + ^ V'-g^ + ^ ^kV^ Qr + 

r=l i=l r=l r=l , , 

n n n n l^O.lo; 

+ ^ ^ v^y>i Vfey- + ^ ^ v'v ft v^g^ - 

r=l i=l r=l i=l 

n n n n 

- E E E E {^hk ^"^v' Pi Ps - 0^.1 v^v' ft ft g,) 



r=l 5=1 i=l j=l 
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Let's compare formula (8.10) for time derivative (p with our previous notations 
(6.1). It's easy to see that (8.10) can be written as 

n / n \ n 

k=i V i=i / fe=i 

Second term in (8.14) is scalar product of vector W and covector ^. Suppose that 

regularity condition is fulfilled (sec definition 6.1). Then we can introduce same 
scalar product into the formula (8.11) for second order partial derivative ilp: 

n / I ^ n n n 

^ = E^^'-^'^ + EE"^^'-^'^+E/3'^-'- (8-15) 

fe=l fe=l r=l fc=l 

Here we used formula (6.3) written in the following form: 

6'^ = p!: + (8.16) 

Combining formulas (8.14) and (8.15) for time derivatives (p and (p, we obtain the 
equality with right hand side free of scalar product (4 | W) : 



^-^^ = EE"'"^^-^^+E/?'^-'- 

fc=l r=l fe=l , , 

(EE — a — + E^^1-^ • 

fc=l \ i=l s=l s=l J 

For the sake of brevity in further calculations we introduce the following notations: 

Vk = Pk-2^2^ ^ 2^ — ^ — . (8.18) 

i=l s=l s=l 

Quantities 771, . . . , ry^ in (8.18) are components of covector field rj. Using notations 
(8.18), we can simplify the above equality (8.17): 



^-^-^^^ = ^Y.^^Pr'-^k + T.Vk-r''. (8.19) 

fe=lr=l k=l 



Deviation function (p is scalar product of vector r and covector p. Using (8.16), we 
can introduce such scalar product into right hand side of (8.19). Then we get 

^-!^^-^^ = ±±arp^^.^, + ±±,^pr.r''. (8.20) 

fc=l r=l k=l r=l 

Now let's recall the definition 8.1. Suppose that the following equalities are fulfilled: 

n n 

E"'"^' = 0, ^VrPL=0. (8.21) 
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Then deviation function (p satisfies second order ordinary differential equation (8.8) 
witli coefficients A and B determined by vector field ot and covector field r}: 

^=M^^ S=M^. (8.22) 

This means that equations (8.21) are sufficient for weak normality condition for- 
mulated in definition 8.1 to be fulfilled. 

Conversely, suppose that weak normality condition for Newtonian dynamical 
system (7.2) is fulfilled. Combining (8.8) with (8.20), we derive 

n n n n 

A^ + B^ = Y,Y.'^''Pr<k+Y.Y.'^rPl-T\ (8.23) 

/c=lr=l fe=lr=l 

where 

A = A-'^, S = B-MZ). (8.24) 

Substituting (8.14) and (1.13) for ip and (f into (8.23), we obtain 

n n n n 

A{$\W)+C{p\T)=Y,Y.^'Pr-^k+Y.T.flrPk-r''. (8.25) 

fe=l r=l fe=l r=l 

Here for more convenience we introduced the following auxiliary notations: 



^ r=l 1=1 r=l 



Vr = Vr 



Aij2^rV'pi + Qr]. (8.27) 



Now conceptual point ! Looking at formula (8.25), let's remember that for a fixed 
trajectory^ = p{t) of Newtonian dynamical system (7.2) functions T^{t), . . . , r"(t) 
and ^i{t), . . . , ^n{t) satisfy system of first order linear homogeneous ordinary dif- 
ferential equations (see (8.4) and (8.5)). Solutions of these differential equations 
form 2n-dimensional linear space T, while components of vcx;tor-function T{t) and 
covector- function ^{t) for any fixed time instant t = to can be treated as coordinates 
in 1. In other words, linear 1 is isomorphic to direct sum 

1^Tp{M)(BT;{M), (8.28) 

where p = p{to). In regular case (sec definition 6.1) denominator = f7(p, p) 
defined by formula (6.2) is nonzero for p ^ 0. Hence vector W = W(p, p) is 
transversal to null-space of momentum covector p: 

<P = ImP = {X e Tp(M) : (p I X) = 0}. (8.29) 

Therefore tangent space Tp{P) in (8.28) breaks into direct sum 

Tp(M) = (W) e (8.30) 
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where (W) is linear span of vector W. We can write similar expansion for T*{M): 

t;{M)^{p)®QS. (8.31) 
Here (p) is linear span of covector p and W is null-space of W: 

W = ImP* = {y e T;(M) : (y | W) = O} . (8.32) 
Substituting (8.30) and (8.31) into (8.28) we get the following expansion: 

1^ (W)©(p)©£©^ = 2:i©T2. (8.33) 

Expansion (8.33) of the space T generates conjugate expansion of dual space T*: 

T* ^ II © 1*2. (8.34) 

The above equality (8.25) is an equality of functions. However, if wc treat pair of r 
and ^ as a point of 1, then for any fixed t = to both sides of the equality (8.25) can 
be treated as linear functional in T, i.e. they are elements of 1*. Note that left 
hand side of (8.25) is in TJ, while right hand side is an element of T2. Subspaces T* 
and I2 in direct sum (8.34) have zero intersection. Therefore both sides of (8.25) 
are zero (as elements of 1* and as functions of t as well, since t = to bs an arbitrary 
fixed instant of time). Thus we have the equality 

^(^|W)+C(p|t)=0. (8.35) 

Due to expansion 11 = (W)* © (p)* from the equality (8.35) we derive 

^ = 0, C = 0. (8.36) 

Applying (8.26), (8.27), and (8.24) to (8.36), we get B = 0, fik = Vk and come back 
to (8.22). Substituting ^ = and B = into (8.23), we obtain the equality 

n n n n 

k—1 r—1 k—1 r—1 

which breaks into two separate parts equivalent to (8.21). Thus, assuming that 
weak normality condition is fulfilled, we have derived again the equalities (8.21). 
When written explicitly, they form a system of partial differential equations relating 
extended vector field V, extended covector field Q, and extended afiine connection 
r. These equations are called weak normality equations. Just above we have 
proved the following theorem for them. 

Theorem 8.1. Weak normality condition for Newtonian dynamical system, (7.2) 
is equivalent to weak normality equations (8.21) that should be fulfilled at all points 
Q = (P; p) of cotangent bundle T*M, where p 7^ 0. 

Using (8.12), (8.13), (8.18), one can easily write weak normality equations (8.21) 
explicitly. However, we shall not do it, since they are rather huge. 
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9. Additional normality condition. 

Let's proceed with studying normal shift phenomenon assuming that weak nor- 
mality condition for Newtonian dynamical system (7.2) is fulfilled. For this pur- 
pose let's consider some hypersurface a and apply initial data (1.14) to differential 
equations (7.2). As a result we get (n — l)-parametric family of trajectories of 
Newtonian dynamical system (7.2). Let's fix some point po G ct- If ■, ■ ■ ■ 
are local coordinates on u in some neighborhood of the point po and if x^, ... , a;" 
are local coordinates in M in n-dimensional neighborhood of this point, then our 
(n — l)-parametric family of trajectories is represented by functions 

' =x^{i,y^,... ,y"-^), f v\ =Pi(i,2/\--- ,2/""^), 

(9.1) 

^ =a;"(t,yi,... ,2/"-i), [ = Pn(i, 2/\ • • • , Z/""')- 

Comparing (9.1) with (8.1) and looking at (8.2) and (8.3), we see that now we can 
define vector-functions Ti, . . . , t„_i and covector-functions ^i, . . . , (,n-\- Their 
components in local chart are given by the following derivatives: 

dx^ 

r/ = ^, U = ^..Vs. (9.2) 

Vector-functions ti, ... , r„_i with components (9.2) define n — 1 deviation func- 
tions yi, . . . , according to the formula (1.13): 

n 

^i = (Pi(i,y\... ,y"-i) = (p|Ti) = ^T/p,. (9.3) 

s=l 

We assumed that weak normality condition is fulfilled (see definition 8.1). Therefore 
each function in (9.3) satisfies differential equation of the form (8.8): 

(p,=A(t,y\... ,y--^)^, + B{t,y\... ,y^-^)Vi. (9.4) 

According to definition 1.4, in order to have normal shift we should provide van- 
ishing of all deviation functions (pi, . . . , (pn-i in (9.3). Due to the equation (9.4) 
it is sufficient to provide the following initial data for them: 

iPi\ =0, ipA =0. (9.5) 

\t=o lt=o 

First part of initial conditions (9.5) is fulfilled due to initial data (1.14). Second 
part of these initial conditions should be fulfilled at the expense of proper choice 
of scalar function u = u{p) = v{y^ , . . . ,y") in (1.14). In order to calculate time 
derivative 0, we can use formula (8.14). Here it is written as follows: 

n / n \ n 

<^i = l] l]v,y'p, + QJ -T^ + J^w'-U- (9.6) 

8=1 \r=l / s=l 

Then we should substitute t = into (9.6). Vectors ti, ... , t„_i at initial instant of 
time t = form base of coordinate vectors in tangent space to initial hypersurface 
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(7. Initial value of momentum covector p is given by (1-14). Thus we have to 
calculate initial values for covectors 4ii • • • > in formula (9.6): 

^si = ^TiPs = - ris + p -VnTis. (9.7) 
lt=o \t=o <jy 

Substituting (9.7) into (9.6) and then substituting (9.6) into (9.5), we get 

£ = E • V.n. - ^ E E V.^>. + q] ■ rt. (9.8) 

^ s=\ s=l \r=\ J 

In two-dimensional case n = dimM = 2 we have only one variable y = y^. Then 
(9.8) turns to ordinary differential equation with respect to function v = z^(y). 
Normalizing c;ondition (1.3) yields initial value problem for this ordinary differential 
equation, which is always solvable (at least, locally). 

In multidimensional case n > 3 we have several variables y^, ... , y"'~^. In this 
case partial differential equations (9.8) form complete system of Pfaff equations 
with respect to function f = v{y^ , . . . , t/"~^). Initial value problem 

v{po) = I/O (9.9) 

at some fixed point po & cr with local coordinates 2/q, ... , 2/o ^^ is typical for Pfaff 
equations (9.8). However, now it is not unconditionally solvable (even locally). 

Definition 9.1. Complete system of Pfaff equations (9.8) is called compatible if 
initial value problem (9.9) for it is locally solvable for all po £ M and for all ^ 0. 

Let's write Pfaff equations (9.8) formally, denoting by tjji their right hand sides: 

0=Vi(i/,2/\... ,y"-^) (9.10) 
Due to (9.10) we can calculate mixed partial derivatives of v in two different ways 



Equating (9.11) and (9.12), we get compatibility condition for (9.10): 

^ij{,y,y\... ,y^-')=^ji{,,,y\... ,y"-'). (9.13) 

Lemma 9.1. Pfaff equations (9.10) are compatible in the sense of definition 9.1 
if and only if for v ^ left and right hands sides of (9.13) are equal to each other 
identically as functions of n independent variables y^, . . . , and u. 

Lemma 9.1 is standard result in the theory of Pfaff equations. Proof of this 
lemma can be found in thesis [7]. 
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Definition 9.2. Wo say that Newtonian dynamical system satisfies additional 
normality condition if PfafF equations for the function in (1.14) derived from 
initial conditions (9.5) are compatible for any hypersurface cr in M. 

For the sake of convenience we introduce extended covector field U with components 

n 

Us=J2^sV^Pr + Qs- (9.14) 
r=l 

Then Pfaff equations (9.8) are written as follows: 

^,.^j:w.V..n.-^j:u..r-. (.15, 

^ s=l s=l 

Now, relying upon definition 9.2, wc shall derive explicit form of compatibility 
equation (9.13). For this purpose let's calculate partial derivatives (9.11) and (9.12) 
in explicit form. For partial derivatives (9.11) we obtain 

\3 n n 



(9.16) 



" ^ s=l r=l 

s—1 r—1 s—1 r—1 

+ E E - M ' E v.. ( ^ ) • v.. n. - 

5=1 r=l S = l ^ ^ 

5 = 1 ^ ^ 5=1 5=1 

In order to transform (9.16) we need to bring about some preliminary calculations. 
For covariant derivative Vx^ rf in (9.16) we have 

rj 5 n n q2 s n n 

V r? = -^+Wr* t^'t" = I Wf" T''r« 

" r=l g=l ^ r=l 9=1 

Taking into account symmetry of connection components F^^ = F*^, we derive 

Vx.t;-V.^.t/=0. (9.17) 
In a similar way by direct calculations we derive the following identities: 

v.(^)^i:v.(2?)...i;.'(2f).e 



(9.18) 



r=l r=l 
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Formulas (9.18) are special cases of general formula applicable to arbitrary extended 
tensor field X. If ^jlZjl are components of X in local chart, then we have 

n n 
r=l r=l 

And finally, there is an identity for commutator of two covariant derivatives: 



n n n 

Pq 



(9.19) 



[V., v.,]n. = - E E E V < -7 + 

q=l a=l 7—1 

n n n n n n 

+ E E E V ^"^"^ - ^ V ^ "^'^ < • 

r=l a=l 7=1 r=l a=l 7=1 

Combining (9.15) and (9.7), for the quantities ^rj in (9.18) we derive 

n n 
^rj = Y.-Pr Vx.n. - E ^ (9-20) 

Here are components of projection operator P introduced in (6.3). Further 

Its is easy to see that right hand side of this formula simplifies when we introduce 
components of projector operator P. Indeed, we have 

^^TT =-^+E^n^- (9-21) 



Combining (9.20) and (9.21), we obtain the following equality: 

r=l ^ ^ fc=l r=l g=l 

fc=l fc=l r=l q=l 



(9.22) 



In similar way for the first term in right hand side of first equality (9.18) we derive 

r=l r=l q=l 

Now we are able to proceed with transforming the equality (9.16). Note that 
terms symmetric in indices i and j make no contribution to ultimate compatibility 
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equation (9.13). Therefore further we shaU omit them replacing by dots. Taking 
into account (9.17), (9.18), (9.19), (9.22), and (9.23), for (9.13) we derive 



n n n n 



M'EEEE ^^^^^^^^^(^'-'V..n,)(P^V.,n,)- 

fc=l q=l r=l 5=1 



71 n n 



q=\ r=l 5=1 \ m=l 
m=l k=l / q=l r=l s=l \ 



m=l m=l fc=l y 



X 



n n / n 



r=l 5=1 \ m=l m=l 

n n ^ 2^ I ^« + 

m—1 ra—1 / 

r=l 5=1 fc=l g=l \ m=l / 



In deriving the above equality we also used the following quite obvious formulas: 

\ O ) O 02 



0/0 02 02 ' 

'"^^ (9.24) 



o y o ^ 02 

m=l 



Using (9.24), for two summands in right hand side of second equahty (9.18) we get 
E^ (o-j ■^'-^=-^EEE o^ Pr^.,n,+ 

r=l ^ ' r=lfe=lm=l 

+ -E E ^ V,n. - ^ ^ ^^4^ U, + (9.25) 



r=l fc=l fe=l m=l 

k=l k=l q=l m=l 

n n 



E^Hti) -^^ =E-^^. -EE — — (9-26) 



r=l ^ r=l r=l m=l 



Right hand sides of the equalities (9.25) and (9.26) are reflected in the above formula 
for 6ij — Oji. In deriving this formula for Oij — 6ji we have made the following 
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transformations for second term in right hand side of commutator identity (9.19): 

n n n n n n 

E E E V -I = E E E V = 

r— 1 a— 1 7—1 k=l m=l r=l 

n n 71 n V * ' 

= E E E E (^^- (i^.^ V..n,) rj= - ^ C/. < rj=) . 



k=l q=l m=l r=l 



Terms from right hand side of (9.27) are also reflected in the above formula for 
6ij — 6ji. Looking at this formula, we see that compatibility equation (9.13) pro- 
viding compatibility of Pfaff equations (9.8) has the following structure: 



71 n n n 



k=l g=l r=l 5=1 

+ E E E § V..n,) r; - i. f E E E § X (9-28) 

q=l r=l s=l q=l r=l s=l 

n n ^ ^ r 

X (P,« V.^.n,) + . 5: 5: r[ r; = 0. 

r=l 5=1 

Here A""*, BJ, and Crs ai"© components of three extended tensor fields A, B, and 
C respectively. They are given by explicit formulas 

^rs ^ v''|^% (9.29) 

n 77, 

Bl = V'-C/, + ^ ^ 1^^= p„ D^/ - 



m=l fe=l 



- V«W^ + 2^ ^ VsVn 



(9.30) 



m=l 



n 

771—1 



n n / n 



(9.31) 



k=l q=l \ m=l 



Further study of compatibility equations (9.28) with coefficients (9.29), (9.30), 
(9.31) require some information concerning geometry of hypersurfaces in non-metric 
geometry of manifolds equipped with symmetric extended connection F and gen- 
eralized Lcgcndrc transformation (2.2) given by extended vector field V in p-re- 
presentation (see (2.4)). 

10. Geometry of hypersurfaces. 

Let's fix some arbitrary point go = iPo^ p) of cotangent bundle T*M. This means 
that we fix some point po G M and some covector p G T*^ {M). Assume that p / 0. 
Then null-space of covector p is a hyperplane in tangent space Tp^ (M) (see (8.29)). 
Let a be some smooth hypersurface passing through the point po and tangent to 
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null-space of fixed momentum covector p at that point. If n = n{p) is smooth 
normal covector field on a, then for p = we have 

p = !/o ■ n(po), where i^o 0. (10-1) 

Taking constant vq ^ from (10.1), we can expand it up to a smooth nonzero 

function ly = v{p) on hypersurface a (or at least in some neighborhood of marked 
point po on a). Function v{p) satisfies the equality 

v{po) = VQ, (10.2) 

which is just the same as normalizing condition in (9.9). So, we can substitute z/(p) 
into (1.14) and use it for defining shift of a along trajectories of Newtonian dynam- 
ical system (7.2). Now we assume that dynamical system (7.2) satisfies additional 
normality condition (sec definition 9.2). This means that for any hypersurface a 
passing through our marked point po € M and for any nonzero constant in 
(10.2) Pfaff equations (9.8) are compatible. Hence compatibility equations (9.28) 
arc fulfilled. Note that in (9.28) we have explicit entries oi v = v{j)) and implicit 
entries of v through p = j/ • n in arguments of extended tensor fields A, B, and C. 
Moreover, we have implicit entries of v in covariant derivatives 



Vx^-^ = S-EErtr..< (10.3) 

^ k=l r=l 

due to connection components F^^ that depend on momentum covector p = i/ • n 
(see (4.14)). But in (9.28) we have no derivatives of function v = I'ip)- Therefore, 
if we write (9.28) only at our fixed point p = po, wc can replace all entries of z/ 
by normalizing constant from (10.2). Using (10.1), we can express n = n{po) 
through our fixed momentum covector p at the point p = po- 

n = n(po) = — . (10.4) 

In this form n{po) is not too specific property of hypersurface a. It determines only 
tangent hyperplane to a at fixed point p = pq. The only parameters in (9.28) that 
depend on fine structure of hypersurface a at the point po are covariant derivatives 
(10.3). Using them, in [20] we have defined a map /: Tp(cr) T*{M). Indeed, if 
T is some arbitrary vector tangent to a, then t = ■ Ti + . . . + a""^ • t„_i. Let 



/(r) = V^n = ^ ^ {a^ V,,n,) • dx\ (10.5) 

r=l j=l 

It is easy to see that (10.5) defines linear map from tangent hyperplane Tp{a) at 
the point p € a to cotangent space T*{M). We consider composite map 

b = -P*o/oP. (10.6) 



Projection operator P* in (10.6) is a conjugate operator for projector P with compo- 
nents (6.3). Remember that P projects onto the subspace Tp{a) G Tp{M). There- 
fore linear map b: Tp{M) — > Tp{M) is correctly defined by formula (10.6). 
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Linear map b defined by formula (10.6) is associated with second fundamental 
form of hypersurface a. Indeed, let's define bilinear form 

6(X,Y) = (b(Y)|X). (10.7) 

Due to the presence of projection operators P and P* in (10.6) we have 

b{X, Y) = KP(X), Y) = 6(X, P(Y)). (10.8) 

Theorem 10.1. Bilinear form (10.7) defined by linear map (10.6) is symmetric. 

When restricted to tangent space Tp{a) of hypersurface a bilinear form (10.7) 
yields second fundamental form of a. Its components 

Pij =b{n,Tj) (10.9) 

define tensor field in inner geometry of hypersurface a. From (10.8) we derive the 
equality 6(X,Y) = fe(P(X), P(Y)). It means that bilinear form (10.7) and linear 
map (10.6) in outer space are completely determined by components of second 
fundamental form (10.9). Further we need the following theorem. 

Theorem 10.2. Let qo = {po,p) be some fixed point of cotangent bundle T*M 

with p 7^ and let projector P be the value of projector-valued extended tensor field 
(6.3) at this point. Then any symmetric quadratic form b in Tp^(M) satisfying the 
equality (10.8) can be determined by some hypersurface a passing through the point 
Po and tangent to null-space of covector p at this point. 

Theorems 10.1 and 10.2 arc proved in paper [20]. Though these theorems arc 
very important for further study of compatibility equations (9.28), we shall not 
repeat their proofs in present paper. 



11. Additional normality equations. 

As in previous section, let's fix some point pq G M and some covector p ^ 
at this point. This means that we fix some point qo = (po)P) of cotangent bundle 
T*M. Let's fix some arbitrary nonzero constant vq ^ and then use formula (10.4) 
for to define another nonzero covector n 7^ at our fixed point pq. Further, let's 
consider various hypersurfaces passing through the point po tangent to null-space 
of covector n. For each such hypersurface a covector n is normal covector at the 
point po- It can be expanded up to a smooth normal covector field n = n{p) (at 
least in some neighborhood of marked point po). Therefore we can build a into a 
framework of shift construction defined by of Newtonian dynamical (7.2). If this 
dynamical system satisfies additional normality condition (sec definition 9.2), then 
we can choose smooth function 1/ = v{p) normalized by the condition (10.2) and 
such that compatibility equations (9.28) are fulfilled. 

Let Ti, . . . , T„_i be basic tangent vectors of a. Applying linear map (10.6) to 
them, we get a set of n — 1 covectors ©i, . . . , dn-i- In other words, we denote 



Oi = bin). 



(11.1) 
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Let's calculate components of covectors (11.1) in local chart. Using formulas (10.5) 
and (10.6) defining linear map b: Tp{M) T*{M), we derive 

n n 

^" = - E Vx,ng = J2 brs rt- (11.2) 

g=l s=l 

Here hra are components of bilinear form (10.7) in local chart. Comparing (9.28) 
and (11.2), we see that compatibility equations (9.28) can be written as follows: 



(-0)^ E E E E i^su r^) - 

k—1 q—1 r—1 s—1 

n n n j-y,^ n n n j-yj. 

(-of E E E § + (-0)^ E E E § X (11-3) 

q—l r—1 s—1 9—1 r—1 s—1 

X ( 6., r/ ) + .0 E E ^^^^^ r[ r; = . 

r=l s=l 



Varying hypersurface a, we can vary components of bilinear form b in (11.3). In 
particular, we can 1) change the sign of b; 2) get zero quadratic form for b. Due 
to these facts (they follow from theorem 10.2) compatibility equations (11.3) split 
into three separate parts. Now they are written as follows: 

n n 

^^(A--A-)^H^s,=0, (11.4) 

n n n 

E E E(^' V - B'r V) < r/ = 0, (11.5) 

r=l s=l q=l 
n n 

EE(<^"-^-)<^/ = o- (11-6) 



Here in (11.4) we used (11.2) again. Vectors Ti, ... , t„_i depend on the choice 
of local chart on hypersurface a. For a fixed point p = po they can be treated as 
arbitrary n — 1 vectors forming base in tangent hyperplane Tp{a). Similarly, due 
to theorem 10.2 covectors 9i, ... , 9n-i can be treated as arbitrary n — 1 covectors 
in null-space of vector W (see (8.32)). Therefore (11.4) and (11.6) reduce to 

n n 

^^(A--A-)p;pi =0, (11.7) 

r—1 s—1 

n n 

Y,Y.i^rs-Csr)P[P,'=0. (11.8) 
r=l s=l 



In this form equations (11.7) and (11.8) do not depend on any particular hypersur- 
face we used to derive them. 

As for (11.5), we should bring it to similar form independent on a. For this 
purpose remember that we can treat ri, ... , t„_i as arbitrary n — 1 vectors in 
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tangent hyperplane Tp{a). Therefore equation (11.5) reduces to the following form: 



71 n n 



E E E(^' V - V) Pi Pj = 0. (11.9) 

r=l s—1 q—1 

In the next step we rewrite (11.9) in coordinate- free form. It looks Hke 

6(P o B o P(X), Y) = 6(X, P c B o P(Y)) (11.10) 

Here B is linear operator in Tp^ (M) determined by components of extended tensor 
field (9.30). Due to theorem 10.2 the above equality (11.10) means that composite 
operator P » B » P is symmetric with respect to any symmetric bilinear form b in 
Tp(,(M) for which (10.8) is fulfilled. Hence composite operator P°BoP can differ 
from projector P only by some scalar factor A: 

PoBoP = A-P. (11-11) 

This fact is proved in paper [20] . Scalar factor A is given by trace formula 

tr(PoBoP) tr(BoP) " " 



n 71 

n-l n-1 -^-^ n-1 ^ ^ 

r=l s=l 



Formulas (11.11) and (11.12) written in local chart yield required equations 

±±PrB:P!=±±§^P;. (11.13) 



In form (11.13) equations (11.5) do not depend on any particular hypersurface used 
in deriving these equations. 

Equations (11.7), (11.8), and (11.13) taken together form a system of additional 
normality equations. Due to the above notations (6.1), (6.2), (6.3), (9.14), (9.29), 
(9.30), and (9.31) they arc partial differential equations for components of extended 
vector field V and extended covector field Q that determine Newtonian dynamical 
system (7.2). We have derived them assuming that Newtonian dynamical system 
(7.2) satisfies additional normality condition formulated in definition 9.2. Con- 
versely, if additional normality equations are fulfilled, then compatibility equations 
(9.28) turn to identities. Therefore Pfaff equations (9.8) appear to be compatible. 
Thus, we have proved the following theorem analogous to theorem 8.1. 

Theorem 11.1. Additional normality condition for Newtonian dynamical system 
(7.2) is equivalent to the system of additional normality equations (11.7), (11.8), 
and (11.13) that should be fulfilled at all points q = {p, p) of cotangent bundle T*M, 
where p 7^ 0. 

Now suppose that both weak and additional normality conditions are fulfilled. 
In this case all deviation functions (9.3) satisfy second order ordinary differential 
equation (8.8) and we can provide initial data (9.5) for them by proper choice of 
function v = u{p) on any predefined hypersurface a in M. Then all deviation 



36 



RUSLAN A. SHARIPOV 



functions do vanish, and we have normal shift of hypcrsurface a. So, we see that 
weak and additional normality conditions are complementary to each other, and if 
both are fulfilled, we can arrange normal shift of any predefined hypersurface along 
trajectories of Newtonian dynamical system (7.2). 

Definition 11.1. We say that Newtonian dynamical system satisfies complete 
normality condition if both weak and additional normality conditions for this 

system are fulfilled. 

According to theorems 8.1 and 11.1, weak and additional normality conditions 
are equivalent to weak and additional normality equations for parameters V and 
Q determining Newtonian dynamical system (7.2). Therefore complete normality 
condition is equivalent to complete system of normality equations including (8.21), 
(11.7), (11.8), and (11.13). As we noted just above, complete normality condition 
is sufficient for strong normality condition to be fulfilled (see definition 1.5). We 
shall strengthen this result in the next section. 

12. Equivalence of strong and complete normality conditions. 

Part of the statement declared in the title of this section is already proved. 
Indeed, we know that complete normality condition implies strong normality con- 
dition. Let's prove converse implication. Assuming that Newtonian dynamical 
system (7.2) satisfies strong normality condition, we should prove that it satisfies 
weak and additional normality conditions. 

In the first step let's prove that additional normality condition is fulfilled. For 
this purpose let's take some arbitrary hypersurface a with marked point p = Po 
and smooth normal covector field n = n(p) in some neighborhood of marked point. 
Then let's take some nonzero constant vq ^ Q and let's apply strong normality 
condition, which is fulfilled by assumption (see definition 1.5). As a result wc get 
smooth function v = ^{p) normalized by condition (10.2) and such that it provide 
initial data (1-14) for normal shift of a. Due to normality of shift all deviation 
functions (9.3) arc identically zero. Hence initial conditions (9.5) for them are 
fulfilled. Writing (9.5) explicitly, we find that our function v = v{p) is a solution 
for Pfaff equations (9.8). Now, varying constant vq ^ Q m (10.2), we prove that 
Pfaff equations (9.8) arc compatible (sec definition 9.1). Thus, additional normality 
condition is proved (see definition 9.2). 

In the second step we shall derive weak normality condition assuming that strong 
normality condition is fulfilled. This is a little bit more complicated. For this 
purpose we fix some point go = {po, p) of cotangent bundle T*M with p ^ 0. Initial 
point po and momentum covector p at this point form initial data for Newtonian 
dynamical system (7.2). They define a trajectory p = p{t) passing through initial 
point Pq. Null-space of initial covector p is a hyperplane in tangent space Tp^^^M). 
Let's consider various hypersurfaccs passing through initial point pq tangent to this 
hyperplane and denote by a one of them. If n = n(p) is normal covector of this 
hypersurface a, then at the point p = po we have the equality (10.1) that determine 
normalizing constant uq ^ for (10.2). Applying strong normality condition (see 
definition 1.5), we find smooth function v = i'{p) on a in some neighborhood of 
initial point po that provide initial data (1.14) for normal shift of a. Thus, our fixed 
trajectory p = p{t) passing through initial point po appears to be shift trajectory 
among many others. Due to normality of shift all corresponding deviation functions 
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(9.3) on this trajectory are identically zero. Therefore we can write 

(Pi\ =0. (12.1) 
lt=o 

Note that initial conditions (9.5) arc also fulfilled. Therefore we can use formula 
(8.20) for ip in left hand side of (12.1). Here it is written as follows: 

n n n n 

=J2J2''^Pr■^k^\ +J2J2'^rP^-rh . (12.2) 
lt=0 lt=o ^^^^^1 lt=o 

Vectors xi, ... , r„_i at initial instant of time t = form base in tangent hyperplane 
to initial hypersurfacc a. As we noted in section 11, for fixed point p = po they 
can be treated as arbitrary n — 1 vectors in null-space of momentum covector 
p = vq ■ n{po). For components of covectors ^i, . . . , ^„_i at initial instant of time 
t = we can use formula (9.20). Applying (9.20) and (11.2), we get 

n n n 

l^^^-^il =^uoP^y^^n, = -J2''obrsT:. (12.3) 

fc=l 8=1 S=l 

Combining (12.1), (12.2), and (12.3), we derive the equality 

n n n n n 

- E E E ^0 "'^ Prhs r/ + E E • = 0. (12.4) 

k=l r=l s=l k=l r=l 

Just like in (11.3), by varying hypersurface a and by applying theorem 10.2 we can 
break (12.4) into two separate equalities 



/c=lr=l fc=lr=l 

which are equivalent to weak normality equations (8.21). Applying theorem 8.1, we 
find that weak normality condition is fulfillcxi, i.e. strong normality condition im- 
plies weak normality condition. Ultimately, we have proved the following theorem. 

Theorem 12.1. Strong and complete normality conditions for Newtonian dynam- 
ical system (7.2) are equivalent to each other. 

13. Connection invariance. 

Theorem 12.1 is a basic result in the theory of dynamical Newtonian dynami- 
cal systems admitting normal shift, while definition 1.5 is basic definition of this 
theory. Comparing them we see that strong normality condition formulated in def- 
inition 1.5 is applicable either to general Newtonian dynamical system of the form 
(2.3), and to special one given by the equations (7.2). Theorem 12.1 is formu- 
lated only for Newtonian dynamical system (7.2), which implies presence of some 
symmetric extended connection F in M. Theorem 5.1 gives one way to avoid this 
discrepancy. We can use symmetric affine connection (5.8) canonically associated 
with dynamical system (2.3) and by means of this connection we can rewrite (2.3) 
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in form of (7.2) (see formula (7.1)). However, there is another way. Below we shall 
prove that the whole theory constructed in sections 6-13 is invariant under gauge 
transformations changing one connection for another: 

T^j^T'^j+T^. (13.1) 

Here T^^- are components of some symmetric extended tensor field T of type (1,2). 
Applying gauge transformation (13.1) to dynamical system (7.2), we change covari- 
ant derivatives VtPi in left hand side. In order to keep corresponding connection-free 
equations (2.3) unchanged we should change components of covector Q as follows: 

n n 

Qi^Qr-Y.Y.^isPkV'- (13-2) 

k=l s=l 

Having fixed gauge transformations by formulas (13.1) and (13.2), now we shall 
apply them to all normality equations (8.21), (11.7), (11.8), and (11.13) for to 
prove their invariancc under these transformations. First of all note that vector 
field V in (7.2), vector field W introduced by formula (6.1), scalar field given 
by formula (6.2), and projector field P with components (6.3) are invariant under 
gauge transformations defined by (13.1), (13.2): 

V y^ W W, (13.3) 

As for covector field U in (9.14), here we have the following transformation rule: 

n n 

«7«^C/«+^^p™r™W^«. (13.5) 

q=l m=l 

Applying (13.1) to curvature tensors (8.6) and (8.7), we derive 

n n 
s— 1 m— 1 

n n n 

+ T.12p^ + E i^im TTr - T^m T^) + (13-6) 

s^l m— 1 m— 1 

n n n n 

+ E E ^'"^^ - E E p« ^"^Ti, 

s=l m=l s=l m=l 

Dt^^D^-V^T^. (13.7) 

Now we can apply (13.1), (13.2), (13.3), (13.4), and (13.7) to vector field a with 
components (8.12) used in weak normality equations (8.21): 

n n n / n n n 

^ - E E E TiPrV^ + ^'y'p^ + 

i=l r=l s=l \ r=l s=l i=l 
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n n n n n n n 



+ E E E + E E E E p-^ p^^^l- 

r=l s=l i=l r=l s=l i=l m=l / 

n n n n n n / n 

~T.T.T.T. v'-v'=v^>. - E E E '^'■^^^^i - 

r— 1 s~l m=l r—1 i=l \ 5=1 

n n n ri \ 

m=l 5=1 m=l 5=1 / 

n n n n n n / n 

+ E E E E ^''v'PiPsV^+Y.T.iY. ^ ^^'^'^ + 

r=l 5=1 z=l m=l r=l z=l \ 5=1 

n n \ n n n 

+ E E v^y'PiVV'- - E E E ^'^'^ 



5=1 m=l / r=l 5=1 m=l 



Looking attentively at the above formula, we see that almost all terms in right hand 
side do cancel each other. As a result we get the following transformation rule: 

a^^a^ (13.8) 

Formula (8.13) for components of extended covector field /3 is more complicated 
than formula (8.12). Let's simplify it using notations (9.14): 



r=l r=l r=l ,^3_g. 

n n n n ' 

r=l r=l 5=1 m=l 

Formula (13.9) for is still rather complicated. Therefore we perform some pre- 
liminary calculations. Using formula (13.5), we derive 

n n n n 

VrU, ^ VrUs +J2J2Prr^ ^^T^ + E E ^'^ - 

q—1 m—1 9—1 m—1 

71 71 71 71 n n 

- E - E E E ^rs Tl^.W" +Y,Y.P- ^rk X 

k—1 k—1 q=l m=l fc— 1 u=l 



n 71 71 71 71 71 n 

^'U.+Y.Y.Y.PuTruT^.^'^^Y.Y.Y.Y.PuPr^ 

k=\ u=\ q=l k=l u=l q=l m=l 



(13.10) 



n 71 n 71 



k=l u=l q=l m=l 

n n 7% 

v-t/, ^ v'-[/, + t:, w^+yHp-^ + 

g=l ^^=1 m=l 

n n 

+ EE^^™^^"^'"^''- 

q=l m=l 



(13.11) 
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For covariant derivative Vfel^* we have the foUowing transformation rule: 

n n n 

s=l 8=1 m=l 

Combining formulas (13.11) and (13.2), we derive 



71 n n 



J2 ^"Uk Qr^J2 ^'^fe + E E ^kg + E E E ^ 

r=l r— 1 r=l q=l r—1 q—1 m=l 



n n n n n n 



X V'-T,- WQr + J^T.T.Pm K W« - E E E ^'^-^ >< 
r=l q=l m=l r=l s=l u=l 

n n n n n n n n n 

r=l s=l q=l u=l r=l s=l q=l u=l m=l 

n n n n n 

X pm v'-Tfe™ w t:^sV' 

r—l s—1 q—1 u—1 m—1 

In a similar way, combining formulas (13.12) and (13.5), we derive 

n n n n n n n 

r=l r=l r=l s=l r=l s=l m=l 

n n 71 n n n n 

r=l q=l m=l r=l s=l g=l m=l 

n n n n n 

X Ti t:!:^w^v^ +Y,Y.Y.Y.Y.p- p-^ v^v^ w . 

r=l 8=1 q=l u=l m=l 

Transformation rule for V^Qr is similar to (13.10): 

n n n n 

VfcQ. ^ VkQr - E E P» '^fe^r" - E E ^rs ' 

S—1 u—1 8=1 U—1 

n 71 n n n n 

-E^'^rQ^+EEJlp^WsV^ + J:Y.p^ X 

q—1 s—1 u—1 q—1 771—1 q—1 

n n n n n n n 

X V«g. -Y.T.Y.PmTr, Tr's V^-Y.Y.T.T.P-PmX 

S—1 q—1 m—1 s—1 q—1 m—1 u—1 



(13.13) 



n n 71 n 



X Tfe- VT" - E E E T^PuPm T^, Vs v«^^ 

s=l q=l m=l u=l 

For term with curvature tensor in (13.9) we use transformation rule (13.7): 

n n n n n 7% 

E E E QsW^ ^y^HHp-^ - 



r=l 8=1 m=l r=l s=l m=l 
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n n n n n n n n 



r=l s=l m=l r=l s=l q=l u=l m=l 



n n n n n 



+ E E E E E ^'^"^ ^"^^'^ 

r=l s=l g=l u=l m=l 

And finally, for term with another curvature tensor -R^^^ in formula (13.9) we use 
transformation rule (13.6), which is more complicated: 

71 'II n n n '11 

r=l 5=1 m=l r=l s=l m=l 

n n n n n n 

+ E E E ^'"^ - E E E V'^^- - 



r=l s=l m=l r=l 5=1 m=l 

n n n n n 



- E E E E E p« (^r ^^'^ - i^;^ + 

r=l s— 1 m— 1 u— 1 q'=l 



n n n n 



E E E E ( ^^'^ ^kr - T,- Tl W^V^) + 

r—1 5—1 m— 1 g— 1 



n n n n n 



+ EE E EEp"P™ (r« v-'t,™ -t,", v-'t-) w^'-y^. 

T-=l 5=1 m=l g=l u=l 

Combining (13.10), (13.11), (13.12), (13.13) and other above formulas, for compo- 
nents of covector field /3 we derive the following very simple transformation rule: 

n n 
m=l q=l 

In deriving (13.14) we used simplified version of formula (8.12) for vector field ct: 

n n n 

r=l r=l r=l 

n n n n 

r=l r=l 5=1 g=l 

In a similar way, using notations (9.14), we can simpliiy formula (8.18): 

Vk=Pk-2_^ ^ • (13.15) 

5 = 1 

Substituting (13.14) into (13.15) and using (13.4), (13.5), and (13.8), we obtain 



n n n 



% ^ % + E E E^™ ^fc" (13-16) 

5=1 m=l q=l 
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Formula (13.16) moans that in general extended covector field 77 is not invariant 
under gauge transformations (13.1), (13.2). However, if first equation in (8.21) is 
fulfilled, then formula (13.16) simplifies. It turns to 

Theorem 13.1. Weak normality equations (8.21) for Newtonian dynamical system 
(7.2) are invariant under gauge transformations (13.1), (13.2). 

Now let's proceed with our calculations for additional normality equations (11-7), 
(11.8), (11.13). Scalar field A given by formula (9.29) does not depend on Q and 
on connection components T^j . Therefore we have 

-^rs ^ ^rs' (13.17) 

For Bl, using (13.5), (13.7), and (13.11), from (9.30) we derive 



n n n 



B: ^Bl-Y,Y.Y.Pk (V™ - V'-W-™) P« . (13.18) 

m=l q=l k=l 

Now let's substitute (13.18) into (11.13) and use (9.29). As a result we get 

n n n n 

^^^^^^^ 



n n n n n 



r=l s=l m=l q=l k=l 

Looking at (13.19), we see that in general left hand side of the equation (11.13) 
is not invariant under gauge transformations (13.1), (13.2). However, if equations 
(11.7) are fulfilled, then formula (13.19) reduces to 

n n n n 

E E ^« ^ E E Pr b: pi (13.20) 

r=l s=l r=l s=l 

For the expression in right hand side of (11.13) under the same assumption we have 

ee!4^^-^ee!4^^- ^''-''^ 

r=l 5=1 r=l s=l 

This follows from trace formula (11.12) for scalar factor A in (11.11). 

Now let's apply gauge transformation (13.1), (13.2) to tensor field C with com- 
ponents (9.31). By rather huge, but direct calculations we find: 



71 n n 



a. ^ a. + . . . + E E E^« t;, p^ v«c/, - ■ 

fc=l g=l «=1 \ 



n n n 



k=l u=l m=l 
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n n n n n / tutq 

E E E E E ^ru [Pt ^'W- Pv Tl^ - V'^H^- ^ V. Tl^ , . 

k=l 0=1 u=l m=l v=l ^ 



By dots we denoted terms symmetric with respect to indices r and s. They do not 
affect ultimate equation (11.8), therefore we need not keep them in exphcit form. 
In right hand side of the above formula we see three distinct terms with sums. For 
the beginning let's transform second term using formula (8.16): 



EE E vuv>^ n = EE E E^'«^-n'= >< 

k=l u=l m=l k=l u=l m=l q=l 



X 



UsPm + J2Yl E E^*? fi p^] ^ 



w=l fc=l \ ^=1 m=l 



p„ t;, ^ = E E E E^^« ^"^ n jy^^ u,pm. 

k=l u=l m=l q=l 

Substituting this result into formula for Crs and taking into account (9.30), we get 

n n n n n n n n 

a«-a. + ...+EEEp-^Hcn'^.' + EEEEEp«^"feX 

k—l <?— 1 u—1 k—1 <?— 1 u—1 m— 1 v—1 

n 71 n 71 n T/f/'? 

X p,'^ V P. - E E E E E v'^i^™ p. r;,. 

/c=l g=l u=l m=l v=l 

Now we apply the same trick with formula (8.16) to last term in the above formula: 

n 71 n 71 n 



E E E E E^« ^"w- t;. 



k—l 9—1 u—1 m—1 v — 1 

71 n 71 n n 7i 



= E E E E E E^« v'^^- ^ t;, + (13.22) 

k—l q—1 u=l m=l ^=1 a=l 



+ EEEE^^^^%^ EE^^^v'^^'"^''" 

/c=l g=l 11=1 v=l \ a=l m=l / 



P_vTljW^ 



Last term in (13.22) is symmetric with respect to indices r and s. Therefore we 
can omit it adding to those denoted by dots in formula for Crs- 



n n n 



n n n n n n 



(13.23) 



■ E E E E E E^'" T^rk vw- p^p, T-„. 

k=l q=l u=l m=l v=l a=l 
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Now we can find transformation rule for left hand side of the equation (11.8). Using 
formula (13.23) and taking into account equations (11.13) and (11.7), we derive 

n n n n 

E(^" - Pi Pj ^ E E(^- - Pi (13.24) 

r=l 5=1 r=l s=l 

Now, summarizing formulas (13.17), (13.20), (13.21), and (13.24), we see that the 
following theorem is proved. 

Theorem 13.2. Additional normality equations (11.7), (11.8), (11.13) for New- 
tonian dynamical system (7.2) are invariant under gauge transformations given by 
formulas (13.1), (13.2). 

Theorems 13.1 and 13.2 mean that we can apply all normality equations (8.21), 
(11.7), (11.8), (11.13) to Newtonian dynamical system written in connection-free 
form (2.3). For this purpose one should set Qi — in them, and one should 
choose identically zero connection components F^^- = 0, thus replacing covariant 
derivatives V and V by corresponding partial derivatives: 

dx^ ' ^ dpi 

However, in this form normality equations arc not obviously coordinate covariant. 
Each term in them loose transparent tensorial interpretation provided by covariant 
derivatives. Generally speaking, we have a problem of constructing present theory 
in coordinatc-frcc and connection-free form. This is separate problem, it will be 
studied in future papers. 

14. Basic example. 

Let H = H{x^, ... , x", pi. . . . ,Pn) be Hamilton function for some Hamiltonian 
dynamical system. It is given by the following well-known differential equations: 

Using some symmetric extended afHne connection F, we can write Hamilton equa- 
tions (14.1) in terms of covariant derivatives V and V: 

x' = V'H, VtPi = -ViH. (14.2) 

Using the same Hamilton function H = H{x^, . . . ,a;", pi, ■ ■ ■ ,Pn) as in (14.2), we 
define so called modified Hamiltonian dynamical system: 

.• = ^^, V.p. = -,^. (14.3) 

^Ps^'H ^Ps^'H 

s=l s=l 
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Surely wo should choose H such that denominator in (14.3) is non-zero. Dynamical 
system (14.3) is an example (is special case) of Newtonian dynamical system (7.2) 
with vector field V and covector field Q given by formulas 

r = ^^. 0. = -^^. (14.4, 

s=l s=l 

Substituting (14.4) into (6.1) and (6.2), we find 

W' = -V\ n = -1. (14.5) 

Theorem 14.1. Modified Hamiltonian dynamical system (14.3) is an example of 
Newtonian dynamical system (7.2) admitting normal shift of hypersurfaces in the 
sense of definition 1.5. 

We can write modified dynamical system (14.3) in connection-free form: 



dH 
dpi 



^ dH 



Pi 



dH 
dx^ 



V- dH 
t^i ^P^ 



(14.6) 



Theorem 14.2. Modified Hamiltonian dynamical system (14.6) is an example of 
Newtonian dynamical system (2.3) admitting normal shift of hypersurfaces in the 
sense of definition 1.5. 

Differential equations (14.3) and (14.6) are different representations of the same 
dynamical system. Therefore theorems 14.1 and 14.2 formulate the same result. 
This result is not new. It was obtained in paper [19]. In present paper it is built 
into framework of more general theory and forms basic example for this theory. It 
proves that class of dynamical systems considered in this theory is not empty, and, 
moreover, this class is sufficiently large. 

15. Dedication. 

This paper is dedicated to my mother F. M. Sharipova, who taught Mathematics 
for many years in School no. 18 of Karakul, Bukhara region, Uzbekistan. She was 
best in recognizing various constellations on the sky and knew great many of them. 
In dark south nights in Summer, when I was 10 or even younger, she often told 
me about stars, comets, planets, and other thing. Possibly this was why later on I 
have chosen Natural Sciences and Mathematics for my profession. Bright image of 
my mother is ever kept in my memory. 
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